Trapping colloids near chemical stripes via critical Casimir forces 
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We study theoretically and experimentally the solvent-mediated critical Casimir force acting on colloidal 
particles immersed in a binary liquid mixture of water and 2,6-lutidine and close to substrates which are 
chemically patterned with periodically alternating stripes of antagonistic adsorption preferences. These pat- 
terns are experimentally realized via microcontact printing. Upon approaching the critical demixing point 
of the solvent, normal and lateral critical Casimir forces generate laterally confining effective potentials for 
the colloids. We analyze in detail the rich behavior of the spherical colloids close to such substrates. For 
all patterned substrates we investigated, our measurements of these effective potentials agree with the corre- 
sponding theoretical predictions. Since both the directions and the strengths of the critical Casimir forces can 
be tuned by minute temperature changes, this provides a new mechanism for controlling colloids as model 
systems, opening encouraging perspectives for applications. 
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1. Introduction 

The structures and the dynamics of colloidal suspensions are to a large extent determined 
by the effective interactions among the mesoscopic colloids which are mediated by the 
solvent surrounding them. A rather intriguing example of such solvation forces is the 
critical Casimir effect 1 1] which occurs close to a critical point of the solvent. As a soft- 
matter analogue of the quantum-electrodynamic Casimir effect 0,0], the confinement of 
the concentration fluctuations within a near-critical binary liquid mixture by the surfaces 
of the colloidal particles results in an attractive or repulsive effective interaction 
Upon approaching the bulk critical demixing point of the mixture at temperature T = 
T c the critical fluctuations become long-ranged so that the resulting effective interaction 
among colloids can extend well beyond the molecular scale of the solvent; it may act even 
at distances comparable to the size of the colloids. The strength and the spatial range of the 
critical Casimir force depend sensitively on temperature because the correlation length d; , 
which characterizes the extension of the spatial correlations of the critical fluctuations and 
therefore sets the interaction range, diverges as E, ~ ^o\(T — T c )/T c \~ v upon approaching 
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the critical point, T — > T c , where V ~ 0.63 is a standard bulk critical exponent and B,q 
is a solvent-specific microscopic length. The critical Casimir effect is universal in the 
sense that all quantities related to it, such as the various components of the force and the 
associated potentials, can be suitably expressed in terms of scaling functions which are 
identical for all those fluids undergoing the same kind of transition, irrespective of their 
molecular details and material properties, which actually enter only via two independent 
non-universal amplitudes, such as £o> an d the value of T c . 

Besides some indirect quantitative experimental determinations of the critical Casimir 
force acting within wetting films 17|-|9|], its direct effect on a spherical colloid close to a 
substrate and immersed in a critical water-lutidine mixture was measured only recently 
ifToh . These experiments clearly demonstrated that, upon approaching the critical point, 
the critical Casimir force can modify the effective interaction between the colloid and 
the substrate by several multiples of k^T . Moreover, depending on the preference of the 
confining surfaces for lutidine or water — denoted as (+) and (— ) effective boundary con- 
ditions (BC) for the order parameter, respectively — the critical Casimir force between 
the colloid and the substrate is attractive [for (±,±) boundary conditions] or repulsive 
[for (+,— ) boundary conditions] 11011 . In addition, these critical Casimir forces depend 
strongly on the geometrical shape of the confining surfaces so that topographically struc- 
tured surfaces experience also a component of the force which is parallel to the walls flUl . 
whereas nonspherical particles are subject to a torque II 1211 . 

One can combine these features of the critical Casimir effect by using surfaces which 
are chemically patterned on the micron-scale in order to obtain also lateral critical 
Casimir forces 11311 . The critical Casimir interaction between a colloid and a substrate 
which is chemically patterned with stripes of alternating adsorption preferences results 
in a laterally confining potential for the colloid, which can be controlled and reversibly 
switched on and off by varying the temperature of the near-critical solvent I14l4la1 . Re- 
cently, such a substrate has been realized by removing via a focused ion beam a regular 
pattern of equally spaced stripes from the hydrophobic HMDS -coating of a glass sur- 
face [14]. The effective interaction between this chemically patterned substrate and a col- 
loidal particle has then been measured by monitoring via digital video microscopy the 
behaviour of a dilute solution of colloidal particles as the demixing critical point of the 
binary mixture of water and lutidine is approached [14]. In comparing these data with 
the corresponding theoretical predictions 01511 it turned out that the critical Casimir poten- 
tial resulting from a chemical pattern depends rather strongly on its geometrical details 
and it actually provides a sensitive tool to probe these microscopic features, which might 
not be easily accessible otherwise lfl5ll . Indeed agreement between the theoretical predic- 
tions of Ref. lf]~3h and the experimental data of Ref. fl4ll for the spatially varying effective 
potential was found only by assuming that the chemical steps between two subsequent 
stripes within the pattern were not microscopically sharp due to fluctuations inherent to 
the preparation process of the structures 11150 . However, in spite of this improvement in 
the comparison between the shapes of the theoretical and experimental potentials, the 
analysis of the temperature dependence of the corresponding potential depths renders an 
amplitude <^o of the correlation length £ (treated as a fitting parameter) which is about 
twice as large H 1 5fl as the one known from independent previous measurements. This per- 
sisting discrepancy, together with the necessity to test independently the aforementioned 
assumption used in Ref. II 1511 . calls for a more detailed experimental study of the critical 
Casimir force acting in the presence of patterned substrates. 

Colloidal suspensions are not only relevant as model systems in soft-matter physics 
but they may find also applications in integrated micro- and nano-devices. Since critical 
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Casimir forces provide a tool to exert active control over the strength and the direction 
of interactions between colloids in such devices, a thorough quantitative understanding of 
the underlying physics is essential. Recently, it was theoretically predicted that for suit- 
able geometrical parameters colloids can levitate in a stable position above chemically 
patterned substrates L16]. Such a "critical Casimir levitation" may even help in overcom- 
ing the problem of stiction which hampers the functioning of currently available micro- 
and nano-electromechanical devices. Therefore, in order to benefit from this wide range of 
possibilities, a theoretical understanding of the phenomenon and the experimental demon- 
stration of its applicability are essential. 

Here, we report on experimental studies of the critical Casimir interaction between a 
colloid and substrates which are chemically structured by microcontact printing of alkane- 
thiols. We shall see below that, in the present context the substrates prepared by this tech- 
nique turn out to exhibit higher chemical and geometrical resolution than those prepared 
via focused ion beam, which was used in previous experiments lfl4ll . A dilute suspension 
of hydrophilic spherical colloids of radius R which impose (— ) boundary conditions to 
a near-critical water-lutidine mixture is exposed to patterned substrates consisting of al- 
ternating and periodically repeating stripes of width L± with boundary conditions (±) so 
that the periodicity is P = L + + L along the ^-direction, whereas the pattern is trans- 
lationally invariant along the orthogonal y-direction. In order to extend previous inves- 
tigations HE! to a wider range of chemical patterns, we consider here sequences of 
stripes with various periodicities P and stripe widths L . In particular this allows us to 
investigate the case in which the effects of neighboring chemical steps interfere, which 
strongly affects the resulting forces acting on the colloid. As in Ref. lfl4ll we measure 
the laterally varying equilibrium spatial distribution of the colloids by digital video mi- 
croscopy and define an effective potential on the basis of the equilibrium number density 
of the colloids projected along the vertical z-direction (orthogonal to x and y). We analyze 
theoretically this potential and the ensuing surface-to-surface distance z of the colloid 
from the substrate and compare our findings with the experimental data. Our theoretical 
analysis provides accurate information also about the vertical probability distribution of 
the positions of the colloidal particles which, however, cannot be resolved by the kind of 
video microscopy used in the current experiment. For all stripe widths and temperatures 
studied we find quantitative agreement between the measurements and the theoretical pre- 
dictions. Both the shapes and the depths of the measured potentials agree very well with 
the theoretical analysis. Moreover, the uncertainty of the actual position of the subsequent 
chemical steps is significantly reduced by using the present microcontact printing in order 
to structure the substrate. The amplitude <§o of the bulk correlation length, which we esti- 
mate by comparing the depth of the measured potentials with the theoretical predictions, 
turns out to be in quantitative agreement with previous, independent estimates. Thus, our 
detailed comparison of experimental data and theory demonstrates that indeed the full 
benefits of the critical Casimir force can be reaped reliably and may be utilized in poten- 
tial applications, opening a new route for using colloids as model systems or in micro- 
and nano-devices. 



2. Theory 

2.1. Critical Casimir potential 

According to renormalization group theory, in the vicinity of the critical point at T = T c 
the normal and lateral critical Casimir forces as well as the corresponding potential can 
be described by universal scaling functions (see, e.g., Refs. 0, S 17]). These scaling 
functions are universal in the sense that they do not depend on molecular details but only 
on the gross features of the system and of the confining surfaces. These gross features, 
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among which there are the spatial dimension d, the symmetries of the system, and the 



geometry of the boundaries, determine the so-called bulk |4j, |5Q and surface [4, 18-20] 
universality classes of the associated critical point. The relevant thermodynamic proper- 
ties which emerge upon approaching the point of the continuous phase transition can be 
understood and analyzed in terms of the fluctuations of the so-called order parameter 
of the phase transition. For the consolute point of phase segregation of a binary liquid 
mixture <j> is given by the difference between the local and the mean concentration of one 
of the two components of the mixture; thus, it is a scalar quantity and the bulk univer- 
sality class is of the so-called Ising type (see, e.g., Refs. ifiol EH]). For classical binary 
liquid mixtures, the surfaces generically exhibit preferential adsorption for one of the two 
species of the mixture, which results in a local enhancement of the order parameter <j) 
close to the surface. This enhancement is effectively described by symmetry-breaking 
surface fields H, 18-2(J and it is usually denoted by (+) and (— ) boundary conditions 



(BC) corresponding to having a preference for (j> > and <\> < 0, respectively, at the sur- 
face. For the water-lutidine mixture we are interested in, one conventionally indicates the 
preferential adsorption for lutidine and water as (+) and (— ) BC, respectively (see also 
Refs. HIOL |21Q ). In the experimental setup described below, the lower critical demixing 
point of the water-lutidine mixture is always approached from the homogeneous (mixed) 
phase at temperatures T < T c upon increasing T at fixed critical composition of the mix- 
ture. From the experimental point of view it is rather difficult to quantify the strength of 
the adsorption preference exhibited by the different portions of the surfaces. Therefore 
the comparison with theoretical predictions requires assumptions, which can be verified a 
posteriori. In the present case we have qualitative experimental evidence that the chemical 
treatment of the surfaces results in rather pronounced adsorption preferences and therefore 
we shall assume that the surfaces are characterized by the so-called strong critical adsorp- 



tion fixed point 1181 1191. We consider here neither the opposite case of weak adsorption 
at the surfaces (see, e.g., Refs. 122112311 ) nor the effects due to off-critical compositions of 
the mixture, which might lead to a bridging transition analogous to capillary condensation 
(see, e.g., Refs. Jill HQ). 

For the film geometry — in which the binary liquid mixture is confined between two 
parallel, planar, and macroscopicaily extended walls at a distance L — the normal critical 
Casimir force per unit area is given by ifnll 

f^ ±) (L,T) = k B T^k { ^ ±) (L/^, (1) 

where the index (+, ±) denotes the combination of BC (+) and (— ) at the two walls. (In 
the absence of a symmetry-breaking bulk field within the film, which would correspond to 
an off-critical concentration of the mixture, one has ft--) = /(+.+)•) The universal scaling 
functions &(+,±) depend only on the scaled variable L/t,, i.e., on the film thickness L in 
units of the bulk correlation length £ . The functions ±) have been calculated exactly 
in d = 2 [30], for d < 4 using perturbative field-theoretical methods 1311] or effective 



theories 13211 . and in d = 3 numerically via Monte Carlo simulations B3I435I1. At present, 
a quantitatively reliable theoretical determination of both &(+,+) and &(+,-) in d = 3 is 
provided only by Monte Carlo simulations. For a chemically structured wall opposing 
another planar wall forming a film geometry, the critical Casimir force has been obtained 
within mean-field theory (MFT) [13] and from Monte Carlo simulations 



For a spherical colloid with (— ) BC opposite to a chemically patterned substrate with 
alternating (+) and (— ) BC the critical Casimir potential has been calculated both nu- 
merically within full mean-field theory (corresponding to d = 4) and by resorting to the 
so-called Derjaguin approximation exploiting the full knowledge of ±) in the film ge- 
ometry for d = 3 [16]. Accordingly, the critical Casimir force and the critical Casimir 
potential for a colloid close to a chemically patterned wall can be expressed in terms of 
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universal scaling functions which take the geometry into account and which depend on 
the following scaled quantities: 

• = z/B,, where z is the surface-to-surface distance between the colloid and the sub- 
strate, 

• A = z/R, where R is the radius of the spherical colloid, 

• S = x/^/Rz, where x is the lateral coordinate of the center of the colloid such that x = 
corresponds to the colloid being located opposite to the center of a stripe with (— ) BC, 

• II = P I \[Rz„ where P is the periodicity of the stripe pattern along the ^-direction, 

• A = L-/P, where L is the width of the stripes with (— ) BC. 

In particular in d = 3 the critical Casimir potential <I>c can be written as fl6ll 

<P C (L ,P,x,z,R,T) = £ B r-tf (A,n,S,0,A), (2) 

z 

where # is a universal scaling function. The realization of the present geometrical setup 
of a sphere facing a plane still represents a challenge for lattice based Monte Carlo simu- 
lations. Accordingly, up to now it is not possible to obtain accurate numerical data for this 
three-dimensional geometry and therefore one has to rely on the Derjaguin approxima- 
tion (DA), based on the assumption of additivity, in order to calculate approximately the 
critical Casimir potential (see Ref. lfl6ll ). Within the Derjaguin approximation, the scaling 
function # [Eq. (f2])] can be expressed in terms of the known scaling functions £(±,-) for 
the film geometry 1 1611 : 

«(A,n, 8 ,e,A-> o) - «■+-><«> + *-■(«» + V-,(«> - «,-,(») a(l ^ @h 

(3) 

where 

/CO 
d/}(/j-i)/r 3 * (±) _)03e) (4) 

are the scaling functions of the critical Casimir potential of a colloid in front of a homo- 
geneous wall as calculated within the Derjaguin approximation ifioll . and 

G)(A,n,S,0) = l + £ {co s (S + n(« + |),0)~co s (S + n( ? i -§),©)) (5) 

n=— °° 

with co s given by 

-7i ds — (TT&im — M (® 1+^/2 ) 

«<™w=*± C>w -«,-,(») ' <6) 

where M(0) = kr + _\(®) —kr__\(&). For = and 0^1 there are analytic expres- 
sions for (O s (see Eqs. (6) and (7) in Ref. lfl5ll . respectively). The accuracy of the Derjaguin 
approximation has been checked numerically within mean-field theory and it has turned 
out that this approximation describes quantitatively the actual behavior of those numer- 
ical data which correspond to > 4 for A < 1 and n > 0.5, as well as to < < 4 
with A < 0.3 and n > 2 11611 . In the experiments discussed further below the correspond- 
ing values of 0, A, and n vary within these ranges so that, assuming that the previous 
quantitative conclusions extend to d = 3, we expect the Derjaguin approximation to be 
quantitatively reliable. 
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2.2. Background forces 

In addition to the critical Casimir force due to the critical fluctuations of the solvent, the 
colloidal particles of the suspension are subjected to additional effective forces which are 
characterized by a smooth and rather mild dependence on temperature. Typical back- 
ground forces acting within the colloidal suspensions of present interest are due to 
(screened) electrostatic and dispersion interactions and to the gravitational field. In a first 
approximation, which neglects possible mutual influences of these forces, their total po- 
tential is given by the sum of the corresponding contributions: (i) the electrostatic potential 
<J> e i, (ii) the gravitational potential <J> g , (iii) the van der Waals interaction <I > v dw> and (iv) 
the effective critical Casimir potential <J>c- 

Electrostatics — This interaction originates from the fact that, due to the formation 
of charge double-layers, the surface of the colloids and of the substrate acquire a sur- 
face charge once immersed in the liquid solvent. As a result, the polystyrene colloids of 
radius R = 1.2/im and (— ) BC immersed in water-lutidine mixtures experience an elec- 
trostatic repulsion from the substrate. The screened electrostatic potential of the colloid 
at a surface-to-surface distance z from a homogeneous substrate with (±) BC is well ap- 
proximated by 



*A,±(z)/(k B T) =exp{-K(z-z±)}, (7) 

where fc 1 is the screening length and Zq describes the strength of the electrostatic repul- 
sion from the substrate with (±) BC. Although the values of these parameters are deter- 
mined by the surface charge of the colloid, the dielectric constant of the mixture etc. [38], 
here they will be treated as fitting parameters of the actual experimental data for <I> c i ±. In 
experimental conditions similar to the present ones as far as the mixture and the colloids 
are concerned, one finds /c" 1 ~ 12nm [HI ED, El andzg ~ 0.1 -O^m Jlol fTlEllEil 



as typical values at T ~ T c ~ 307 K. In view of a possible difference between Zq and Zq , 
the resulting electrostatic potential of a colloid close to a patterned substrate such as the 
one described above depends on the lateral position x of the colloid. Such a dependence 
can be accounted for theoretically within the Derjaguin approximation (by repeating the 
calculations described in Appendices A.2, B, and C in Ref. jlol for Oc), which is ex- 
pected to be particularly accurate due to the exponential decay of <& e l ± in Eq. Q) as a 
function of the distance z from the substrate. For this electrostatic potential & e i(x,z) one 
therefore finds (for a colloid facing the center of a (— ) stripe at x = 0) 



<S>el,+(z)+4>el,-(z) , ^el+(z)-^cL-(z) ^, v ... 

$> c \{x,z) = 1 (8) 

where 

oo 

Q(x) = l+ £ {erf([x + «P-L_/2]/A)-erf([x + «P + L_/2]/A)}, (9) 



with A = V2Rk~^ — 0.17/im. Q. and <I> c i depend on the geometric parameters L and P 
describing the pattern. 

Dispersion forces — For the particular choice of materials and conditions used in the 
present experiments^van der Waals forces turn out to be negligible compared with the 
other contributions 1 10 , 15 , EO, HH] . In addition, the temperature dependence of the dielec- 
tric permittivity of the water-lutidine mixture is negligible (see footnote 4 in Ref. LI 5]) and 
therefore no significant changes are expected to occur in Ovdw an d ^ci within the range 
of temperatures explored experimentally. 

Gravity — Due to buoyancy the colloid immersed in the solvent above the patterned 
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substrate experiences a gravitational potential given by 

4k i 

*g(z) = (PPS - PWL)g Y R Z = GZ ' (10) 

where pps — 1.055g/cm 3 and p\yL — 0.988g/cm 3 [40] are the mass densities of the 
polystyrene colloid and of the water-lutidine mixture at the critical composition and 
near T c , respectively, and g ~ 9.81m/s 2 is the gravitational acceleration. Accordingly, at 
T ~ 307K one has G ~ 1 . HksT / and therefore, compared to the other contributions, 
it turns out that the gravitational potential depends rather mildly on the distance z because 
it varies only over a few k B T on the relevant length scale of a few microns. The expres- 
sion in Eq. ( fTOl assumes that the colloidal particle of mass density ppg is floating in a 
homogeneous medium of mass density Pwl- However, the laterally varying adsorption 
preferences of the substrate induce the formation of alternating water-rich or lutidine-rich 
regions close to the surface of the patterned substrate, which laterally alter the resulting 
mass density of the solvent as a consequence of water and lutidine having different mass 
densities. This implies that the effective gravitational constant G acquires a dependence 
on x. In addition, the preferential adsorption of the colloid, with the ensuing formation of 
an adsorption profile around it, can lead to a modification of the effective density ppg of 
the colloid itself. However, on the basis of our estimates, we expect all these effects to be 
negligible for the present experimental conditions 11511 . 

2.3. Total potential 

The total potential <I> of the sum of the forces acting on the colloid is given by 

&(x,z,T) = 3>c(L-,P,x,z,R,T)+3> c i(x,z)+<5> g (z), (11) 

where the theoretical expressions for the individual contributions are given by Eqs. (O, 
([8]), and (fTOl . Here and in the following we do not indicate the dependence of <I> on L , 
P, and R, because the values of these parameters are fixed for each individual experiment. 
Figure Q] shows the total potential &(x,z,T) of a single colloid with (— ) BC opposite to 
a chemically patterned substrate, as a function of both x and z and for three values of the 
temperature T close to the critical value T c . These three values correspond to different 
correlation lengths as indicated in the figure. The gray area in the x-z plane indicates 
the vertical projection of the stripe with (— ) BC, the center of which corresponds to x = 0. 
In Fig.[TJthe white part of the x-z plane corresponds to the projection of the stripe with (+) 
BC, the center of which is located at x = P/2 = 0.9/^m. The potential of the forces acting 
on the colloid is translationally invariant along the ^-direction which is not shown in Fig.Q] 
<&(x,z,T) in Fig. Q] has been calculated by using geometrical and interaction parameters 
which mimic the actual experimental conditions and by using values of the correlation 
length E, which are experimentally available. As anticipated above, Fig. Q] clearly shows 
that the gravitational tail of the potential, which characterizes <E>(;t,z, 7") at large values 
of z, is indeed rather flat on the scales of IcrT and of hundreds of nm. As a consequence 
of thermal fluctuations — which cause the colloid to explore a region of space within 
which the total potential <I> differs from its minimum by a few ksT — the particle is 
expected to display large fluctuations Az ~ kgT /G along the vertical direction. At small 
particle-substrate distances electrostatic forces are responsible for the strong repulsion of 
the colloid from the substrate. In Fig. Q] we allowed for a lateral inhomogeneity of the 
electrostatic potential (i.e., Zq ^ Zq), which might occur due to different surface charges 
on the different stripes. This is clearly visible in panel (a) of Fig. [Q which corresponds 
to a rather small value of the correlation length ^ so that, within the range z — 0.1/im 
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z Iflm] z [urn] z [fim] 

(a) t, = 5nm (b) t, = 22nm (c) t, = 26nm 



Figure 1. Total potential &(x,z,T) of a colloid with radius R = 1.2^m opposite to a chemically patterned substrate with 
L_ = 0.9/Jm and P = 1.8/lm for three temperatures corresponding to (a) t, = 5nm, (b) ^ = 22nm, and (c) | = 26nm, 
respectively. The electrostatic potential corresponds to Zq = 0.12tlm, Zq = 0.08/im, and If -1 = 12nm (see the main text for 
details). At small separations z. < zj, the colloid is strongly repelled from the substrate due to electrostatics. The effective 
gravitational potential associated with buoyancy is rather flat, with a spatial slope G ~ l.\2ksT / flm. Accordingly, the 
surface-to-surface particle-substrate distance z exhibits large thermal fluctuations of the order of kg T/G ~ 1 tlm as long as 
t, is small. Upon approaching the critical point [from (a) to (c)], S, increases and a deep, local potential minimum arises 
rapidly as a function of t, above that part of the substrate (indicated by the shaded area) with the same preferential adsorption 
as the colloid. The colloid is eventually confined in this deep potential well at a distance z — Zq with almost no vertical 
fluctuations. The black lines correspond to cuts through the potential surface at constant values of z and x, respectively. 

Fig. (Ha) refers to, the contributions of the critical Casimir force are negligible. In this 
case, the x-dependent electrostatic contribution dominates at small values of z, whereas 
the laterally homogeneous gravitational potential dominates at larger distances. However, 
upon approaching the critical point [Figs, [lib) andQJc)], the correlation length B, increases 
and the critical Casimir force acting on the colloid builds up; it is repulsive within the 
region with (+) BC whereas it is attractive within the region with (— ) BC. In the latter 
case the behavior of the colloid is eventually determined by the electrostatic repulsion 
and the attractive critical Casimir force, whereas in the former case it is determined by 
gravitation and the repulsive Casimir effect. Above a certain threshold value of ^ , which 
depends on the specific choice of the various geometrical and physical parameters, a very 
deep and steep potential well develops rapidly close to the stripe with (— ) BC, which 
therefore confines the vertical motion of the colloid at much smaller values of z than 
before with very limited thermal fluctuations of the particle-wall distance. In contrast, the 
vertical repulsive critical Casimir force, which the colloid experiences above the stripe 
with (+) BC, pushes it further away from the surface, but the corresponding fluctuations 
of the vertical position z (still limited only by the gravitational tail) are not significantly 
affected. Thus, the full theoretical analysis of the various forces at play reveals a rather 
interesting energy landscape which is strongly temperature dependent. 



2.4. Measured potential 

In order to measure experimentally the total potential <I> of the forces acting on a colloid 
a very effective approach consists in monitoring the equilibrium Brownian motion of a 
single particle and inferring from the sampled probability distribution function P(x,y,z) x 
exp(— <&(x,z,T) /(ksT)) the potential as <&/(1cbT) = — \nP + const, where x and y are 
the lateral coordinates of the projection of the colloid center onto the substrate surface. 
This approach forms the basis of total internal reflection microscopy, which has been 
used in Refs. iQH EH l23ll to study critical Casimir forces. Alternatively, one can study 



a colloidal suspension which is sufficiently dilute so that the inter-particle interaction is 
negligible. In this case, the mean equilibrium number density p(x,y,z) of the colloids at 
position (x,y,z) is proportional to the single-colloid probability distribution function P 
and therefore it is again given by p(x,y,z) 06 exp(— <&(x,z,T) /ksT). In the experimental 
setup described below, the positions of the centers of the colloids are monitored via a 
digital video camera positioned below the substrate. Accordingly, the surface-to-surface 
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distances z of the colloids from the substrate are not resolved and the camera records 
only the projected number density pp(x,y) = Jq dz p(x,y,z). Due to the translational 
invariance of the chemical pattern along the y direction of length / 3> P,R, the density 
pp(x,y) can be conveniently projected further onto the x-axis, resulting in an effective 
number density p(x) = /~ ! J 'dv pp(x,y) = l^ 1 j^dy J °°dz p(x,y,z), which depends only 
on x. This projection increases the statistics and therefore the accuracy with which this 
projected density can be determined experimentally. Subsequently, an effective potential 
V(x) (up to an irrelevant additive constant) can be associated with p(x) such that p(x) <x= 
exp(— V {x) / {ksT)). (Note that due to the thermal fluctuations of the colloids along the 
vertical direction, even if one knows the average distance z avg (x) of the colloid from the 
substrate at a certain lateral position x, the effective potential V (x) is not simply given by 
<J>(x,z av g(x), T) +const, as it was implicitly assumed in Ref. uM-) The measured potential 

8V(x) = V(x) - V(P/2) = -k B T\n{p(x)/p{P/2)) (12) 

is eventually defined such that it vanishes for a colloid opposite to the center of a repulsive 
(+) stripe at x = P/2. We emphasize that 8V (x) contains universal ingredients stemming 
from the scaling function associated with <I>c (see Eq. (fTTT )) as well as nonuniversal con- 
tributions due to <I> c i and <!>„. 



2.5. Non-ideal stripe patterns 

Due to the preparation process (see below) the actual position x = x s (y) of each chemical 
step separating two adjacent stripes might vary smoothly along the y-axis. This variation 
affects the measured effective potential V of the colloids as long as it occurs on a length 
scale which is comparable to or smaller than the typical distance £ msc j along the v-axis 
which each particle explores during the acquisition of the images by the camera. The im- 
ages aquired during the experiments allow one to estimate such a mean-square displace- 
ment £ mS( j to be of the order of tens of jxm lfl4[[l5 . 41]. The projection along the v-axis, 



which yields the density p(x), effectively causes a broadening of p(x) compared to the 
case of straight (ideal) chemical steps with x s (y) = const. In addition, locally a smooth in- 
trinsic chemical gradient of the step leads to such an effect, too. For illustration purposes, 
we first consider a single chemical step, which is ideally located at x = and which gen- 
erates a potential <J>(x,z,J). (This reasoning can be extended to the periodic chemical 
pattern we are presently interested in by assuming additivity of the forces, see Ref. Illotl 
for details.) In order to estimate the effect of these variations of the actual position of the 
step along the v-axis, we assume that the local position x s (y) of the step does not change 
significantly along the v-axis on the scale of the radius R of the colloid, so that on this 
scale it can be considered as ideal and therefore generates a potential <I>(x — x s (y),z,T). 
For x s (y) we assume an effective Gaussian distribution p(x s ) along the y-axis, with zero 
average and standard deviation Ax. Accordingly, the projection / _1 J ' dy along the y-axis 
turns into j x oo dx s p{x s ) and affects the resulting projected density p(x) and the resulting 
potential V(x) [see Eq. (fT2ll. 



2.6. Particle-substrate distance 

The theoretical knowledge of the total potential <fr(x,z, T) [see Eq. (fTTIl l allows one to 
predict the particle-substrate distance z as a function of the lateral variable x and tempera- 
ture, a quantity which is not accessible to the experiments discussed below. As anticipated 
above, the rather small value of G in Eq. (fTOb is responsible for rather large fluctuations 
of the particle-substrate distance z around the position z = Zmin(-x) at which the potential 
<J>(x,z,r) has a minimum as a function of z for a fixed lateral position x of the particle 
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and which corresponds to the position of mechanical equilibrium. In the presence of these 
large fluctuations it is convenient to consider the x-dependent mean particle-substrate dis- 
tance z av g{x), which is determined by the probability distribution function P of the colloid, 
i.e., by the total potential a^| 

z**g(x) = - r -[ dzzexp{-<t>(x,z,T)/(k B T)}, (13) 
N(x) Jo 

where N(x) = J$dz exp{— <I>(x,z,r)/(fcgr)} is the normalization. In order to describe 
the thermal fluctuations of the vertical position of the colloid it is convenient to consider 
the probability P<(z;x) that for a fixed lateral position x the colloid has a surface-to-surface 
distance from the substrate smaller than a given z: 

P<(z;x) = J- f dz' exp{-<t>(x,z',T)/(k B T)}. (14) 
N(x) JQ 

In order to generalize the notion of "standard deviation" to the present case of an asym- 
metric distribution of the particle-substrate distances at fixed lateral position, we define 
a lower value z\ow{x) and an upper value z upp (jc) of the particle-substrate distances such 
that 



P<(zkm(x);x) =0.159 and P<(z upp (x);x) = 1 -0.159, (15) 

so that the probability of the colloid to be at a distance z with zi ow (x) < z < z U pp{x) 
is ~ 68%, whereas the probability to find it at distances smaller (larger) than zi ow (x) 
(z upp (x)) is ~ 16%; these two properties define the standard deviation for a Gaussian 
distribution. 

Figure [2] shows the behavior of the colloid with (— ) BC above a chemically patterned 
substrate for the two lateral positions (a) x = and (b) x = P/2 at which the colloid is 
floating above the center of a (— ) and of a (+) stripe, respectively, as a function of the 
correlation length In Fig.|2]the values of the geometrical parameters (P, R, L_) are cho- 
sen to correspond to the actual experimental conditions, whereas the parameters Zq gov- 
erning the electrostatic repulsion from the substrate are varied within a range which was 
determined by previous independent experiments. Panel (a) clearly demonstrates that the 
fluctuations of z for a colloid opposite to an attractive stripe decrease significantly upon in- 
creasing the correlation length d; above a certain threshold value t, * which depends on Zq ■ 
This reflects the emergence of the deep potential well shown in Fig. Q] which results from 
the competition between an increasingly attractive critical Casimir force and a repulsive 
electrostatic repulsion, the former being always overwhelmed by the latter around z = Zq ■ 
Indeed, for small values of £; the average position Za, V g{x = 0) is typically determined 
by the competition between the electrostatic repulsion and the gravitational part, such 
that Za, V g{x = 0) ~ Zq +ksT /G ~ 1/im with fluctuations of the order of ksT /G ~ 1/j.m. 
On the other hand, for larger values of one has z U pp,avg,iow,min(* = 0) ~ Zq and there 
are only very small thermal fluctuations of the particle-substrate distance, at most a few 
tens of nm. Depending on the relative strength of the electrostatic repulsion and the criti- 
cal Casimir attraction, one can have z m in < ziow if the former dominates the latter, i.e., at 



2 Note that the Derjaguin approximation holds only for distances which are small on the scale of the particle size (a detailed 
analysis of its applicability for the system under consideration is given in Ref. 1 16]). However, in Eq. 1131 also large values 
of z occur. But at these large particle-wall distances the critical Casimir force as well as the electrostatic force are negligibly 
small compared to the gravitational force, so that using this approximation is nonetheless not detrimental. In principle, the 
integration in Eq. U3t is limited by the vertical extension of the experimental sample cell of around 200^im. However, 
due to the gravitational contribution to the potential, de facto no colloidal particle moves out of the vertical field of view 
of the imaging objective. Thus the integration in Eq. )13> can be taken to run up to infinity without quantitatively relevant 
consequences because contributions from large z are strongly suppressed. 
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Figure 2. Distances of the colloid from the substrate, for a fixed lateral position (a) x = and (b) x = P/2, as a function 
of the bulk correlation length Here, P = 1.8/xm, R = 1.2/im, L_ = 0.9/im, whereas z upp> z aV g> Zlow> ar, d Zmin indicate 
the upper, average, lower, and potential-minimum distances of the colloid, respectively (see the main text). The solid lines 
correspond to the choice Zq = Zq = 0.09/im, the dashed lines to Zq = 0. \2jlm and Zq = 0.09/lm, and the dotted lines to 
Zq = 0.09/im and z.q = 0.15/xm. As expected, the fluctuations of the particle-substrate distance z for a colloid opposite to 
an attractive stripe [x = 0, panel (a)] decrease significantly upon increasing the correlation length t, , due to the fact that 
the particle is mainly localized around the deep potential well which forms as a consequence of the action of the critical 
Casimir force. In fact, z U pp, aV g,iow,min — Zq for large The actual value of Zq strongly affects the behavior of the particle, 
both in determining the values of z U pp, aV g,iow,min close to the critical point and in setting the threshold value of £, at 
which one observes such a sharp transition towards strong spatial confinement for t, > £*. From panel (a) one infers, e.g., 
that ~ 17nm and ~ 27nm for the solid and dashed curves, respectively. On the other hand, the dependence on Zq is 
not pronounced and indeed in panel (a) the dotted lines practically coincide with the solid ones. For a colloid opposite to a 
repulsive stripe [jc = P/2, panel (b)], due to the repulsive nature of both the electrostatic and the critical Casimir force and 
the very weak gravitational attraction one has z m i n < zi ow and the average position is typically of the order of 1/lm with 
fluctuations of the same order. As anticipated, also in this case the actual value of Zq is not very relevant, in particular for 
large values of £, . Analogously, the behavior of the colloid at x = P/2 is not affected by the choice of Zq and in panel (b) 
the dashed lines practically coincide with the solid ones. 



small values of E, or the opposite at large values of t, . From Fig.|2a) one concludes that the 
choice of Zq strongly affects the behavior of the particle, both in determining the values of 
£upp,avg,iow,min an d in setting the threshold value E, * of £ above which the particle becomes 
strongly confined, whereas the dependence on Zq is negligible for the behavior at x = 0, 
because Zq controls the electrostatic interaction with the adjacent stripe. Analogously, the 
behavior of a colloid at x = P/2 as shown in Fig. E£b), i.e., opposite to a repulsive stripe 
is not affected by the choice of Zq ■ However, for this configuration, also the actual value 
of Zq does not affect significantly the resulting behavior of the particle-substrate distance 
at x = P/2, in particular for large values of Indeed, due to the repulsive nature of both 
the electrostatic and the critical Casimir force and the weak gravitational attraction, the 
average position is typically of the order of z a vg(* = P/ty ^ Zq +ksT /G ~ JcbT /G (with 
an additional linear contribution oc £ for large values of £ ) allowing fluctuations of the 
order of ksT/G ~ 1/im. 

In Fig. [3] the particle-substrate distance (characterized via z av g,upp,iow) and the position 
of mechanical equilibrium z m i n are reported as functions of the lateral coordinate of the 
colloid within the period P and for two different temperatures, i.e., two different values of 
§ . In Fig. [3l the choice of the parameters corresponds to actually experimentally accessi- 
ble values. Figure [3] clearly shows that for t, = lOnm (dashed lines) the particle-substrate 
distance is laterally constant and the colloid does not react to the presence of the chemical 
pattern on the substrate, apart for a possible effect due to a change in the electrostatic in- 
teraction (which here is taken to be the same on the different stripes). On the other hand, 
for = 20nm (solid line) the colloid is strongly attracted to that part of the substrate with 
the same preferential adsorption (in Fig. [3] indicated by a shaded background), as a con- 
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Figure 3. Particle-substrate distance z (characterized via z U p Pl avg,low> see the main text) and position of mechanical 
equilibrium Zmin as functions of the lateral coordinate of the particle for t, = lOnm (dashed curves) and £, = 20nm (solid 
curves). Here, Zq = Zq = 0.09/jm, P = 1.8,11m, R = 1.2,um, and L_ = 0.9 jlm. For small values of the correlation length 
£, the various characteristic distances are almost independent of the lateral coordinate, whereas above a certain threshold 
value of | (see also Fig. [2} the colloid opposite to the attractive stripe is strongly confined close to the wall at a distance 
Z — Zq with almost no fluctuations. 

sequence of the emerging critical Casimir forces. Indeed for \x\ < L_/2 the particle is 
abruptly localized at a distance z — Zq , which is primarily set by the electrostatic repul- 
sion and which corresponds to the position z m i n (0) of the minimum of the potential, with 
almost no thermal fluctuations. In the region above the repulsive stripes (in Fig. [3] indi- 
cated by a white background) and within the temperature range considered here, instead, 
the repulsive critical Casimir force pushes the colloid only slightly further away from the 
substrate, by a distance of the order of the increasing correlation length E, . But the am- 
plitude of the thermal fluctuations of the particle-substrate distance are barely affected by 
the onset of the repulsive critical Casimir force. 

The analysis of the particle-substrate distance shows that the behavior of the colloidal 
particle and the resulting potential are drastically influenced by the strong attraction of 
the colloid close to an attractive stripe. Accordingly, the electrostatic repulsion from this 
stripe, i.e., the value of Zq , affects significantly the total potential and has therefore to 
be considered carefully in the comparison between theoretical predictions and the exper- 
imental data. On the contrary, the actual value of Zq as well as the actual value of K do 
not significantly affect the resulting behavior of the potentials as long as they are within 
the range appropriate for the present experiment which can be inferred from previous, 
independent measurements 112 ill . 



3. Experimental setup 

In the present context of critical Casimir forces, surfaces with a periodic pattern of alter- 
nating stripes with opposite adsorption preferences for the two components of the binary 
mixture of water and lutidine were obtained by micro-contact printing (;llCP) of alkane- 
thiols. Since there is a vast literature covering this technique 142144411 here we shall de- 
scribe it only briefly. Stamps were prepared by casting poly(dimethylsiloxane) (PMDS) 
onto a master which was topographically structured by a lithographic process. After the 
PMDS is cured one obtains an elastic stamp which exhibits the negative structure of the 
master. In the next step the stamp is wetted with a ImM ethanolic solution of nonpolar 
HS(CH2)i7 — CH3 thiol. When the excess liquid is removed, the stamp is brought into 
mechanical contact for several seconds with a glass substrate coated by a 30nm layer of 
gold. This results in hydrophobic regions the geometry of which corresponds to the pro- 
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trading parts of the stamp. Finally, the substrate is dipped into a ImM ethanolic solution 
of polar HS(CH2)n — OH thiol for several minutes in order to render the remaining bare 
gold surface hydrophilic. With this technique we were able to produce periodic arrays of 
hydrophilic and hydrophobic stripes of widths between 0.9/im and 3/im over a typical 
total extension of 0.5cm 2 . The typical lateral edge resolution of the chemical structures 
is of the order of 50 — lOOnm. In previous studies IU4I1 the chemical patterns were cre- 
ated by removing locally a hydrophobic monolayer from a hydrophilic surfaces with a 
focused ion beam (FIB). Although under ideal conditions the resolution of FIB is within 
the range of several nanometers, in the case of non-conducting surfaces we observed dis- 
tortions/deflections of the ion beam due to an electrostatic charging of the glass surface, 
leading to deviations of the chemical steps from straight lines. This causes considerable 
discrepancies between the measured critical Casimir forces acting on particles near such 
surfaces and the corresponding theoretical predictions, which assume straight and sharp 
boundaries for the stripes [15]. Since such charging effects are basically absent in jliCP, 
here we have an improved control of the geometrical structure of the imprinted chemical 
pattern and therefore much better agreement with theory (see below). 

In order to probe the critical Casimir forces which act on a colloidal particle exposed 
to the patterned substrates prepared as described above, we used polystyrene (PS) probe 
particles with radius R = 1 .2/im. Their surface charge, as provided from the manufacturer, 
is 10/iC/cm 2 which renders them hydrophilic, realizing (— ) boundary conditions. The 
particles were dissolved in a critical water-2,6-lutidine mixture which has a lower critical 
demixing point at a lutidine mass fraction of c£ = 0.286 and a critical temperature of T c = 
307K. Particle positions were monitored by digital video microscopy which allows one 
to track the projection of their centers onto the substrate plane with a spatial resolution of 
about 50 nm. For each measurement the particle trajectories were recorded for 40 minutes 
with a frequency of 4 frames per second. Within the field of view of 290/imx225/^m 
around 20 to 40 particles have been monitored, corresponding to a projected density of one 
particle per square of the size of 25x25 to 17 x 17 particle diameters. Thus, the particle 
density was sufficiently small to exclude the presence of critical Casimir forces among 
neighboring colloids, such that only their interaction with the patterned substrate is probed 
during the measurements. 

The temperature of the sample was controlled by the procedure described in detail in 
Ref. [14]. With this setup, temperatures close to T c could be stabilized within 10 mK 
over several hours. In contrast to temperature changes, which could be resolved within 
mK accuracy, larger errors in the determination of absolute temperatures (in particular the 
measured critical temperature) occur. This is due to the fact, that the metal-resistance ther- 
mometer (PtlOO) could not be placed within the sample cell but it was instead attached 
outside but close to the field of view. Indeed, the thermometer measures a temperature 
r out lower than the actual temperature T inside the sample cell. We associate the temper- 
ature r out = T° ut with the critical temperature T = r c exp of the solvent in the sample at 
that particular temperature for which critical opalescence is observed when shining a laser 
beam into the sample cell. This leads to systematic errors on the absolute temperature of 
T° nt and r c exp of the order of 50 mK. Using the assumption that T — T out = const (this 
constant could not be determined and is inter alia dependent on the ambient temperature), 
we have measured the temperature difference AT = T° ut — T out and identified it with 
r^exp _j = Since the comparison with the theoretical predictions depends crucially 
on the actual value T c ~ 307K of the critical temperature, in the analysis below we account 
for such a possible systematic error by considering AT* = T C — T c exp as an additional fit- 
ting parameter, where T c ~ 307K is the actual critical temperature of the water-lutidine 
mixture. 
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4. Comparison of theory with experiments 

The total potential <I> of the forces [Eq. CCD)] acting on the colloidal particle has been 
calculated on the basis of the Derjaguin approximation with k/± \ [see Eqs. (GQ) and ((6])] 
obtained from Monte Carlo simulations. In the following we shall use for the scaling 
functions fc(± _) of the critical Casimir force between two planar walls with (±, — ) BC 
the numerical estimate referred to as "approximation (/)" in Figs. 9 and 10 of Ref. J34ll . 



We have checked that different [34] or more recent and accurate 113311 estimates for &(±,-) 
actually lead to essentially the same effective potentials, the only difference being a small 
additional overall shift of the resulting fitted value of the critical temperature (see further 
below). A detailed analysis in spatial dimension d = 4 suggests that the Derjaguin approx- 
imation we have used in our theoretical predictions for d = 3 should be rather accurate in 
describing the actual potential of the colloid within the range of parameters experimen- 
tally studied here [16]. As anticipated above, in order to predict the effective potential V 
we need to fix also the value of the parameters which determine the electrostatic inter- 
action. For the comparison between theory and experiment we fix the screening length 
K" -1 to the value K~ l = 12nm which has been reported from independent measurements 
on the same system (see, e.g., Ref. [21]). In order to fit our theoretical predictions to the 
experimental data, we vary instead the unknown values of the parameters Zq [Eqs. (0) 
and ©] within the range 0.08 — 0.15/^m, which can be reasonably expected on the ba- 
sis of previous measurements on homogeneous substrates 12111 . However, the results for 
8v (0) are hardly affected by the particular choice of Zq so that in the following we keep 
it fixed at Zq = 0.09 ;Um. The amplitude f£ of the correlation length has been determined 
by independent experiments as f£ ~ 0.20±0.02nm (see, e.g., Tab. Ill in Ref. 112 ill ) so 
that we vary within the range 0.18 — 0.22nm in order to obtain the best fit to the ex- 
perimental data by the theoretical scaling functions. Moreover, as mentioned above, the 
experimental uncertainty in the determination of the absolute value of the critical tem- 
perature T c exp is taken into account by considering as an additional fitting parameter the 
shift AT* = T C - r c exp (see the end of Sec. of up to |AJ C *| ~ lOOmK. (The values of 
AT* may be different for the individual stripe widths L which characterize the substrates 
investigated in independent experimental runs.) However, the relative uncertainty in the 
determination of the temperature within a single experimental run for a given stripe width 
L is smaller than ± lOmK (see above). Thus, for the comparison carried out below, we are 
left with AT* , Zq, and Ax (see Subsec. |231 ) as fitting parameters which, however, are 
all limited to rather small ranges of values in order to be in accordance with independent 
and previous experimental results. 



4.1. Depth of the potential 

The depth of the measured laterally confining potential is given by AV = —8V(0) > 
[see Eq. (fT2l)1. Figure Ufa) shows AV as a function of the deviation AT = r c exp — T of 
the temperature from the experimentally determined critical demixing temperature r c exp , 
for various stripe widths L and measured in independent experimental runs as described 
above. The periodicities P of the various patterns the stripes belong to are determined from 
the photolithography mask and are given in Tab. [1] The potential depth AV is determined 
as the difference of the potential between its value at the center of a repulsive stripe and 
at the center of an attractive stripe, carrying an uncertainty of ±0.15^7 near the critical 
temperature T c ~ 307K. Upon approaching T c , for a stripe width L < 1.5/im the effects 
of two adjacent chemical steps interfere, which results in an effectively reduced potential 
depth compared to the case of a very wide stripe (see the dashed curves in Fig. ffib); 
in Fig. @Ja) the experimental uncertainty in determining the critical temperature r c cxp is 
responsible for the relative displacements among the various curves). 
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Figure 4. Depth AV of the effective potentials of a colloid close to patterned substrates, with periods P = L- + L + and (— ) stripe widths L_ 
given in Tab. [T] as a function of the temperature deviation from the critical temperature. The symbols represent the experimental data which 
are affected by an experimental uncertainty of ±0.15^7" for the potential depth and of ±().01K for the temperature differences between the 
various data points belonging to one and the same value for L_ . 
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(a) AV as a function of AT = T°* p — T. (For better visibility, the data for 
L_ = 3.25/im are shifted by 0.1K along the AT axis, as indicated by the red 
arrow.) The solid lines represent the best fit of the theory to the experimental 
data for each individual value of L_ with fixed parameters K and Zq which 
do not affect the resulting behavior significantly as discussed in the main 
text. £j" has been varied within the reasonable range ().18nm< < 0.22nm 
known from the literature. The values of all parameters corresponding to the 
solid lines are given in Tab. [T] Clearly, for each individual L_ the actual 
critical temperature T c is shifted by AT* = T c — r r cxp from the value r c cxp 
determined experimentally, with AT* being within the experimental accu- 
racy A7" c * < 0.1K. The dashed curves represent a common fit to all data. 
For this latter fit we have assumed £^ and z,q to be the same for the various 
independent experimental runs, independently of the corresponding value of 
L_ , whereas AT,.* = T c — T c Gxp (besides Ax as determined below) is used for 
adjustment to the data for each individual L_ . The best fit for all these pa- 
rameters renders & = 0.22nm, Zq = 0.11/im, and the values reported in 
Tab.0 




(b) AV as a function of A7" + AT* = T C — T using the val- 
ues of AT* for the best common fit (see dashed lines in (a)) 
given in Tab.|2] This plot takes into account the experimen- 
tal uncertainty of measuring the absolute value of T r cxp . 
The data points, which have been shifted accordingly, ba- 
sically fall on top of each other within their error bars. The 
dashed lines correspond to the dashed ones shown in (a) 
and represent the theoretical predictions. For L_ > 1.60jim 
the various curves are almost indistinguishable from each 
other, whereas for L_ < 1.30iJm the potential depth is ef- 
fectively reduced due to the interference of the effects of 
two neighboring steps and due to the fact that the steps are 
non-ideal. 



The effect of non-ideal chemical steps may result in an effectively reduced potential 
depth as well (see, for example, Fig.HJb) and, c.f., Fig.[6£b)). Accordingly, in comparing 
the experimental data for AV with the theoretical predictions we allow for an uncertainty 
in the local position of the boundary between adjacent stripes, as described in Subsec. l231 
However, as will be shown below, we are able to determine Ax rather precisely from the 
lateral variation of the potential. The values of Ax, which yield the best agreement and 
which are used for the comparison shown in Fig. HJ are given in Tab. [TJ On the basis 
of our theoretical analysis, it turns out that the temperature dependence of the potential 
depth becomes independent of the stripe width when the latter is sufficiently large. Adopt- 
ing for the geometrical and physical parameters the values corresponding to the present 
experimental conditions, this is expected to occur for L > 2/im. From Fig. [4] one can 
see, however, that the experimentally determined data for individual experimental runs 
are shifted along the temperature axis with respect to each other, which reflects the uncer- 
tainty of up to lOOmK in measuring the absolute value of the critical temperature. On the 
other hand, within a single individual experimental run corresponding to a certain stripe 
width, the temperature difference between the various data points can be measured with 
the high accuracy of less than lOmK. Upon comparing the experimental data shown as 
symbols in Fig. 0] with our theoretical predictions we take this into account by introduc- 
ing, for each value of L , a possible shift AT* (L_ ) = T c — T c cxp (L_ ) between the actual 
critical temperature T c and the value J c exp determined in that particular experiment. Ac- 
cordingly, for the data corresponding to a certain L , the actual distance from the critical 
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Table I. Values of the parameters for which best agreement is obtained between theory and the data 
for 8v{x) for each individual experimental run, corresponding to a single value of L_ (see Fig. |4la)). 
The values K" -1 — 12nm and zt — 0.09/im are fixed because their choice does not affect significantly the 
resulting theoretical curves in Fig. [4] 
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a Most of the experimental data for the potential AV discussed here are not influenced by the 
actual value of P as long as P > 2L_ . Accordingly, we assume the periodicity to be the one 
determined by the inscribed photolithographic mask pattern reported here. 

b Width of the stripes of the photolithographic mask pattern. Due to the /J CP stamping pro- 
cess we estimate the uncertainty of the actual width of the thiol stripes to be up to ±0.5/im. 



c Value of the stripe width for which best agreement between theory and experiment is ob- 
tained. 

d Value of Ax for which best agreement between theory and experiment is obtained (see 
Figs.[5]and[6j. 



Table 2. Values of the parameters for which best agreement is obtained between theory and all 
experimental data together (see Fig.l4lb)). so that the values of £ + and z.q are the same, and AT* (in 
addition to Ax determined from the shape of the potentials; see below) is the only parameter allowed 
to vary for the various stripe widths L_ . As in Tab.m the values of k: -1 — 12nm and z.q — 0.09um 
are fixed. (For a description of the parameter values P, L*l xp , and L_ see the footnotes in Tab.fTJ) 
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point is given by T c - T = AT + AT* (L_ ) . 

The dashed curves in Fig. 0] correspond to the common fit to all experimental data 
which leads to the fitting parameters and Zq taking the same values %q = 0.22nm and 
Zq = 0. 1 ljum for all L as obtained from the least-square method. On the other hand, the 
individual temperature shifts AT*(L-) are given in Tab.[2]for each data set corresponding 
to a single value of L . The values of & and Zq are both in agreement with independent 
previous findings [21]. Figure |4jb) shows the depth of the potentials as a function of 
AT + AJ C *(L_) = T C — T, i.e., shifted along the temperature axis by AT*(L_) as given in 
Tab. [2 This accounts for the uncertainty in determining experimentally the absolute value 
of the temperature (see the discussion at the end of Sec. [3]), and indeed, as expected for 
significantly wide stripes L S> £, , the data overlap with each other within the error bars, 
reflecting that in this limit AV is independent of L . 

The solid lines in Fig. HJa) correspond to fits of the theoretical predictions to each 
individual experiment dealing with a specific stripe width. Distinct from the previous 
common fit, here for each stripe width L we allow for a variation of the values of ^ and 
Zq in addition to AT C and Ax. On the basis of a least-square fit, best agreement is obtained 
for the values given in Tab. [TJ which agree with those obtained in independent previous 
studies i2lh within the corresponding experimental accuracy. 



4.2. Shape of the potentials 

Figures[5]and[6]show the total potential of the forces acting on the colloid as a function of 
its lateral position for various temperatures and for two stripe widths. Symbols represent 
the experimental data, whereas the solid and dashed lines are the corresponding theoret- 
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Figure 5. Lateral variation of the effective potential SV(x) (see Eq. 1 1 21 ) for a colloid opposite to a chemically patterned substrate and immersed in the 
water-lutidine binary liquid mixture at its critical concentration for various temperatures r c cxp — AT. Symbols indicate experimental data, whereas the 
lines are the corresponding theoretical predictions for ideal (dashed lines) and non-ideal (solid lines) stripe patterns. The parameters used to calculate the 
theoretical curves are given in Tab.[T] 




(a) L- = 3.25,um. From top to bottom the critical point is approached and 
the measured temperature deviations AT are 0.3, 0.18, 0.15 (0.153), 0.145 
(0. 14), 0. 1 3 (0. 1 28), 0. 1 2, 0. 1 1 (0. 1 1 3), and 0. 1 K with an uncertainty of up 
to ±0.0 IK with respect to each other. If indicated, the values in brackets are 
corrected values (but compatible within the experimental inaccuracy) which 
have been used for evaluating the theoretical predictions. 
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(b) Same as (a) for L_ = 2.25iJm and AT 
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ical predictions with ^ , z , and AT* fixed to the values reported in Tab. [T] which have 
been determined from the fit of the depth AV of the potential. In addition, for a few cases 
and as indicated in the figure captions, we use the leeway provided by the experimental 
uncertainty of ±10mK for the temperature value. The dashed lines in Figs. [5] and [6] are 
based on the assumption that the stripe patterns are ideal whereas the solid ones refer to 
non-ideal patterns for which we fitted the parameters Ax (see Subsec. l231) in order to ob- 
tain the best agreement between theoretical predictions and the set of experimental data 
at distinct temperatures. The resulting values of Ax for the various widths of the pattern 
are reported in Tab. Q] In fact, the broadness of the transition regions across the chemi- 
cal steps of 8V (x) between its extremal values is affected practically exclusively by Ax 
among the fitting parameters. We have checked that the other parameters and potential 
additional effects such as polydispersity and weak adsorption cannot account for the dis- 
crepancies of the slopes of the dashed curves and the experimentally determined ones 
(see also Ref. lfl5ll ). Thus, Ax and L are determined de facto independently of the choice 
of Zq, and AT* and in the following they can be regarded fixed upon variation of 
the values of these latter parameters. In particular, this is important for narrow stripes as 
shown in Fig. [6l where even the depth of the potential depends on L and Ax, in contrast 
to broader stripes. 

As anticipated, in the case of rather large stripe widths L = 3.25/im and L = 2.25/im 
shown in Fig. [5] and for the temperatures considered here the effects of two neighboring 
chemical steps do not interfere and the resulting potential across the center of the attrac- 
tive chemical stripe is flat. Accordingly, from the comparison between the experimental 
data for the shapes of the potentials and the corresponding theoretical prediction we can 
infer the width L of the stripes — now treated as a fitting parameter — independently 
of its initially assumed value L * 13 determined by the width of the photolithography mask. 
Indeed, the actual width of the stripe imprinted by the mask may differ by up to ±0.5jUm 
from L e * v characterizing the mask itself. It is only L e * p which is controlled during the 
preparation process 14111 . On the other hand, the actual value of P is the same as the 
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Figure 6. Effective potential SV(x) (see Eq. )12t ) for various temperatures Tc 
experimental data, whereas the lines are the corresponding theoretical predk 
parameters used to calculate the theoretical curves are given in Tab.[T] 




-2-1012 
x [,um] 

(a) Same as Fig.0a) for L_ = 1.60/xm and AT = 0.3, 0.23 (0.218), 0.21 
(0.198), 0.19 (0.182), 0.17, 0.15 (0.153) K. 



— AT and relatively small widths of the attractive stripes. Symbols indicate 
us for ideal (dashed lines) and non-ideal (solid lines) stripe patterns. The 
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(b) Same as (a) for L_ = 1 .30(im (left panel) and L_ = 0.9/im (right panel) 
with AT = 0.3, 0. 12, 0. 1, 0.09 (0.083), 0.08 K for L_ = 1 .30^m and AT = 
0.3, 0.18, 0.16, 0.145, 0.13 K for L_ = 0.9^m. 



one of the photolithography mask because the periodicity is not affected by the stamping 
process. (Note that, within the parameter range investigated in Fig. [5j the effects of neigh- 
boring chemical steps do not influence each other and therefore the actual value of the 
periodicity P of the pattern does not affect the determination of L_.) In agreement with 
this observation, the comparison with the theoretical predictions leads to fitted values L 
which are within the aforementioned range indicated in Tab. [TJ From Fig. [5j we find that 
the experimental techniques described above lead to rather sharp chemical steps between 
two stripes with an uncertainty of only Ax < 0.15/im. 

The comparison for narrower stripe widths L < 2/xm is shown in Fig. [6] Even for these 
cases, in which the effects of two neighboring chemical steps do interfere, the experi- 
mental data agree very well with the theoretical predictions obtained from the Derjaguin 
approximation. Whereas for L = 1.6/xm and L = 1.3/im we have to take into account 
an uncertainty Ax of the position and associated with the shape of the chemical steps of 
about 0.2/im, for the smallest stripe width L = 0.9/^m the chemical pattern turns out to 
be almost ideal with \Ax < 0.1/im (see Tab.Q]). 



5. Conclusions and Outlook 
5.1. Summary 

We have presented a detailed theoretical and experimental study of the effective solvent- 
mediated forces acting on colloids suspended in a near-critical binary liquid mixture of 
water and 2,6-lutidine and close to chemically patterned substrates. In contrast to ear- 
lier investigations [14] we have obtained a chemical pattern of stripes of different widths 
and with laterally alternating adsorption preferences for either lutidine [(+) boundary 
condition] or water [(— ) boundary condition] by monolayers of thiols deposited on gold- 
coated substrates via microcontact printing. The solute colloidal particles — 2.4/xm di- 
ameter polystyrene spheres — were rendered hydrophilic, corresponding to (— ) boundary 
conditions, due to their surface charge. Based on digital video microscopy our main ex- 
perimental findings are the following. 
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1 . At the fixed critical composition of the water-lutidine mixture, upon raising the tem- 
perature T towards its lower critical demixing point at T = T c , lateral forces acting on 
the colloidal particles arise gradually. These critical Casimir forces cause the colloidal 
particles to be attracted to those chemical stripes which have the same adsorption pref- 
erence as the particles and to be repelled from those stripes with the opposite adsorption 
preference. 

2. These lateral and normal contributions to the total effective forces (compare Sub- 
sec. [2T2l» are negligible at temperatures T more than a few hundred mK away from the 
critical value T c but they increase significantly upon approaching T c . Eventually, they 
lead to laterally confining potentials for the colloidal particles with potential depths of 
several kgT (see Fig.H)). These potentials can be reversibly tuned by minute temperature 
changes. 

3. Across a chemical step the critical Casimir potentials vary rather abruptly between 
two plateau values on a length scale of « 0.8/^m (see Fig. [5]). This indicates that the 
microscopic structures of the chemical steps formed by imprinted layers of alkanethiols 
with different endgroups are much sharper than those created previously by a focused 
ion beam [ 14f| . 

4. For rather narrow stripes of widths < 2/im the effects of two neighboring chemical 
steps interfere and consequently reduce the potential depth and lead to rounded shapes 
of the effective potentials (see Fig. [6]). 

These experimental observations can be consistently interpreted in terms of the occurrence 
of the critical Casimir effect, and it is possible to quantitatively compare the measured 
potentials with the corresponding theoretical predictions. We have derived the effective 
potentials within the Derjaguin approximation (see Subsec. [2Tl and Ref. lfl6ll ) based on 
universal scaling functions for the critical Casimir force between two plates with (±, ±) or 
(+, — ) boundary conditions as obtained from Monte Carlo simulation data. (Differences 
in estimates for these scaling functions as obtained from various Monte Carlo simulations 
do not affect our results significantly.) In contrast to the present experimental measure- 
ments we are able to theoretically analyze the spatially fully resolved critical Casimir 
potential acting on a colloid. The resulting energy landscape for a colloidal particle is 
rather rich due to the interplay of several relevant forces, which we have taken into ac- 
count (see Fig. [T]). Typically, thermal fluctuations lead to large fluctuations of the lateral 
position and the distance between the colloidal particles and the substrate (see Fig. |2]). 
Upon approaching T c , however, strong normal and lateral critical Casimir forces abruptly 
localize the colloids very close to and above stripes exhibiting the same adsorption prefer- 
ence as the particles (see Fig. [3]). In calculating the effectively one-dimensional, projected 
potentials for a single colloid, as they are obtained experimentally, we also take into ac- 
count the possibility that the chemical pattern is not ideal, in addition to other experimen- 
tal uncertainties. Not only the theoretical predictions for the potential depths agree with 
the experimental data (see Fig. 0]) but also the correlation length , as determined from 
our comparison, follows rather well the expected universal power-law behavior and the 
associated non-universal amplitude <§o is in agreement with previous independent experi- 
mental determinations (see Tab.Q]). Moreover, the shapes of the potentials as a function of 
the lateral position of the colloid show good agreement between the theory and the exper- 
iments (see Figs.[5]and[6]>. Since critical Casimir forces probe rather sensitively the details 
of the geometry of the patterns 11511 . we find that the chemical steps obtained by the ex- 
perimental method used here are sharp, with deviations of only a few percent in terms of 
the stripe width. From the detailed comparison with rather narrow chemical stripes (see 
Fig. [3b)) we infer that even for such cases the Derjaguin approximation describes the 
actual behavior quite well, as expected from a corresponding theoretical analysis within 
mean-field theory fl^l . We conclude that the quantitatively successful comparison be- 
tween the experimental data and the theoretical predictions demonstrates the significance 
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and reliability of the critical Casimir effect for colloidal suspensions and reveals a new 
means of using them as model systems in soft-matter physics or in applications in inte- 
grated micro- or nano-devices. 



5.2. Perspectives 

"After many failures, involving many fluid samples, I did succeed in levitating a hard 
platelet." 

— Robert Evans. 

Since quantum-electrodynamic Casimir forces acting on the scale of nm and /im are typ- 
ically always attractive, they are responsible for stiction occurring in micro- and nano- 
mechanical devices. According to a generalization of the Earnshaw theorem, for quan- 
tum Casimir forces acting on conventional materials in vacuum, stable levitation is not 
possible [45]. Therefore, in order to overcome the problem of stiction using the quantum- 
electrodynamic Casimir effect the system must be immersed in a fluid and the bodies have 
to be designed very specifically |4g, 47]. Obviously however, for a system immersed in 



a fluid it is natural to use solvent-mediated forces such as the critical Casimir forces to 
control its behavior. Since one can obtain attractive and repulsive critical Casimir forces 
by suitable surface treatments as discussed here, a properly chosen combination of chem- 
ical substrate patterns may easily lead to levitation of a colloidal particle at a stable height 
above the substrate flfill . One finds that this critical Casimir levitation not only can be 
achieved without involving any other additional force but that the stable equilibrium dis- 
tance between the wall and the colloid can be tuned steeply by varying the temperature 
flfill . Recently JlHl it was predicted that the ranges of geometrical parameters for a pattern 
of laterally alternating stripes (as discussed here), which allow for critical Casimir levita- 
tion of colloidal particles, are rather large. Moreover, critical Casimir levitation is robust 
even in the presence of other forces and of thermal fluctuations fl6ll as they occur typ- 
ically in colloidal suspensions, so that an experimental realization of this effect appears 
to be possible. These results show that chemical patterning of substrates allows one to 
design the critical Casimir effect over a wide range of properties. 

"Ions are the invention of the devil." 
— Robert Evans. 

By adding salt to the solvent of the colloidal suspension the electrostatic repulsion be- 
tween the substrate and a colloid is strongly screened; this provides the possibility to ex- 
plore the spatial variation of critical Casimir forces over a much wider range of distances 
and correlation lengths of the fluid. Moreover, a rather complex behavior of the resulting 
forces acting on the colloid can be expected because electrostatic forces induced by ions 
lead to rather delicate physical features. Recently, preliminary measurements have shown 
a rich and novel behavior of a colloid immersed in a salt-rich solvent upon approaching 



the critical point [48]. The onset of strong attractive forces even several K below the crit- 
ical temperature is observed and, in addition, these forces are found to remain attractive, 
independent of the boundary conditions, (±,±) or (+,—), throughout a large tempera- 
ture range. Currently it is not clear whether these observations can be entirely attributed 
to critical Casimir forces or whether, e.g., the coupling between electrostatic interactions 
and critical fluctuations has to be considered. Thus, using these intriguing effects together 
with a chemical patterning of substrates may open a new field of phenomena with poten- 
tial applications in colloid science. 
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